Search for: All records

Groundwater transit time distributions (TTDs) describe the spectrum of flow‐weighted apparent ages of groundwater from aquifer recharge to discharge. Regional‐scale TTDs in stream baseflow are often estimated from numerical models with limited calibration from groundwater sampling and suggest much younger groundwater discharge than has been observed by discrete age‐dating techniques. We investigate both local and regional‐scale groundwater TTDs in the Upper Middle Loup watershed (5,440 km²) overlying the High Plains Aquifer in the Nebraska Sand Hills, USA. We determined flow‐weighted apparent ages of groundwater discharging through the streambed at 88 discrete points along a 99 km groundwater‐dominated stream segment using³H, noble gases,¹⁴C, and groundwater flux measurements at the point‐scale (<7.6 cm diameter). Points were organized in transects across the stream width (3–10 points per transect) and transects were clustered in five sampling areas (10–610 m in stream length) located at increasing distances along the stream. Groundwater apparent ages ranged from 0 to 8,200 years and the mean groundwater transit time along the 99 km stream is >3,000 years. TTDs from upstream sampling areas were best fit by distributions with a narrow range of apparent ages, but when older groundwater from downstream sampling areas is included, the regional TTD is scale dependent and the distribution is better described by a gamma model (α ≈ 0.4) which accommodates large fractions of millennial‐aged groundwater. Observations indicate: (a) TTDs can exhibit spatial variability within a watershed and (b) watersheds can discharge larger fractions of old groundwater (>1,000 years) than commonly assumed.

We utilized 251 measurements from a recently developed automated seepage meter (ASM) in streambeds in the Nebraska Sand Hills, USA to investigate the small‐scale spatial variability of groundwater seepage flux (q) and the ability of the ASM to estimate mean q at larger scales. Small‐scale spatial variability of q was analyzed in five dense arrays, each covering an area of 13.5–28.0 m²(169 total point measurements). Streambed vertical hydraulic conductivity (K) was also measured. Results provided: (a) high‐resolution contour plots of q and K, (b) anisotropic semi‐variograms demonstrating greater correlation scales of q and K along the stream length than across the stream width, and (c) the number of rows of points (perpendicular to streamflow) needed to represent the groundwater flux of areas up to 28.0 m². The findings suggest that representative streambed measurements are best conducted perpendicular to streamflow to accommodate larger seepage flux heterogeneity in this direction and minimize sampling redundancy. To investigate the ASM's ability to produce accurate mean q at larger scales, seepage meters were deployed in four stream reaches (170–890 m), arranged in three to six transects (three to eight points each) per reach across the channel. In each reach, the mean seepage flux from ASMs was compared to the seepage flux from bromide tracer dilution. Agreement between the two methods indicates the viability of a modest number of seepage meter measurements to determine the overall groundwater flux to the stream and can guide sampling for solutes and environmental tracers.

Stream bank storage effects during floods have received limited attention, despite the significant role of such floods in aquifer water budgets. One reason is the complexity of geometry of the problem, which commonly has been treated numerically. Using a simple model in a domain with moving boundary, a semianalytical solution for bank storage effects is proposed to account for stream stage hydrograph, floodplain slope, and aquifer parameters. The results extend classic solutions by Cooper and Rorabaugh (1963,https://doi.org/10.3133/wsp1536J) for idealized vertical streambanks but applied to realistic floodplain cross sections. The accuracy of the semianalytical solution is verified by a one‐dimensional numerical method and compared to a vertical two‐dimensional variably saturated‐flow numerical model. Comparison indicates that a robust solution is valid for diagnostic analyses of modeling bank storage effects on floodplains. The semianalytical solution is applied to laboratory experiments as well. The results indicate that the present solution provides reasonable estimates of peak timing and head of groundwater flow response in the sloping bank during varying stream stage.